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Abstract. We have analysed SCUBA 850/im images of 
the (near) face-on spiral galaxy NGC 6946, paying par- 
ticular attention to the subtraction of sky signal. A com- 
parison with both 21cm HI and 12 CO(2-l) intensity maps 
reveals a tight correlation between dust thermal emission 
and molecular gas at the kiloparsec level. By means of a 
Monte Carlo radiative transfer model, we convert a B-K 
colour image of NGC 6946 into a map of visual optical 
CN ] depth. The model yields maximum opacities since we as- 
sumc that any increase in B-K colour, with respect to the 
\ disk edge, is attributable solely to extinction by dust. The 
resultant map of visual optical depth relates well to the 
O | distribution of neutral gas (HI+H 2 ) and implies a global 
gas-to-dust ratio of 90 (this value is a lower limit). There 
is no significant radial variation of this ratio: this can be 

■ understood, since the gas content is dominated by far by 
I \ the molecular gas. The latter is estimated through the CO 

■ emission tracer, which is itself dependent on metallicity, 
similarly to dust emission. In the absence of a more objec- 

■ tive tracer, it is not possible to derive the true gas-to-dust 
ratio. By comparing the radial profile of our visual opti- 
cal depth map with that of the SCUBA image, we infer 

j an emissivity (dust absorption coefficient) at 850/im that 

■ is 3 times lower than the value measured by COBE in 
the Milky Way, and 9 times lower than in NGC 891. We 

" o ■ view this very much as a lower estimate, however, given 
our initial assumptions in deriving the visual opacity, and 
the possibility of underestimating the large-scale submm 
| emission, the effect being more severe for the nearly faec- 
es ■ on orientation of NGC 6946. A decomposition of the spi- 
ral structure half way out along the disk of NGC 6946 
suggests an interarm optical depth of between 1 and 2. 
These surprisingly high values represent 40-80% of the vi- 
^ sual opacity that we measure for the arm region. 
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1. Introduction 

Modern submillimeter (submm) and far-infrared (FIR) 
facilities such as the Submillimeter Common User 
Bolometer Array on JCMT (SCUBA), Institut de Ra- 
dioastronomie millimetrique (IRAM) and the Infrared 
Space Observatory (ISO) have undertaken an invalu- 
able role in re-evaluating the radial distributions and 
global masses of galactic dust in nearby spiral galaxies 
(e.g. Guelin et al 1993; Guelin et al 1995; Neininger et 
al 1996; Alton et al 1998a; Siebenmorgcn et al 1999). 
Nevertheless, detailed information about how grains are 
distributed with respect to the spiral arm structure, or 
how they relate at the kiloparsec level to the various gas 
phases (HI, H2) has been less forthcoming (Knapen & 
Beckman 1996; Seigar & James 1998; Rand et al; 1999; 
Alton et al 2001). A major obstacle, in this respect, is 
the limited spatial resolution available with instruments 
such as ISO (~ 100" at 200^m) - barely sufficient to 
resolve the spiral arms in even the closest galaxies. At 
wavelengths closer to 1 mm, ground-based arrays such 
as SCUBA and IRAM, afford considerably better spatial 
resolution (<~ 10" ) but emission arises from the Rayleigh- 
Jeans tail in this regime and, consequently, surface 
brightnesses are very low (usually ~ 10 mJy/beam). The 
sparse emission is particularly problematic for face-on 
spirals and it is precisely galaxies with this kind of 
inclination that we need to study if we are to understand 
how gas and dust relate to the spiral density wave. 



masses inferred from optical extinction/reddening (Block 
et al 1994; Alton et al 1998b). Nevertheless, at the local 
level (~ 1 kpc) it is less clear how grain thermal emission 
relates to optical absorption lanes, star-forming regions, 
the diffuse interstellar medium (ISM) and, of course, the 
various gas phases (HI,H 2 ,HII). 

In a recent paper, Alton et al (2001) collated and 
examined the most up-to-date estimates of spiral disk 
opacity. What became obvious in this study is that 
relatively little is known about how dust is distributed 
with respect to the spiral arms and, in particular, whether 
interarm grain material is prevalent or not. The intcrarm 
region is particularly important because it occupies more 
than 50% of the disk and may therefore be crucial in 
determining the ability of spiral disks to block background 
light (Heisler & Ostriker 1988; Fall & Pei 1993; Gonzalez 
et al 1998; Masci & Webster 1999). Furthermore, the 
distribution of dust with respect to the spiral density 
wave may hold clues as to how gas is compressed along 
the spiral arms to form new stars. We should not forget 
that half of the metals in spiral galaxies are encapsulated 
in grains and consequently galactic dust provides a 
valuable 'fossil record' of previous star-formation within 
the disk environment (Whittet 1992; Williams & Taylor 
1996; Edmunds & Eales 1998). 



At the global level several important trends have 
emerged from recent submm/mm observations of spirals. 
For example, Dunne et al (2000) have shown that the 
total 850/zm flux density correlates more tightly with 
integrated CO(1-0) line emission than with the 21cm HI 
line emission. The mid-IR emission (7 to 15 microns) 
detected with ISOCAM does not correlate well with HI 
and CO either (e.g. Genzel & Cesarsky, 2000, Frick et al 
2001). Total dust masses that include submm/mm data 
now begin to show signs that they are comparable to dust 



Table 1. Basic properties of NGC 6946. 



Parameter 
Type 

R25 radius 

R.A.(J2000) 

Dec. (J2000) 

Distance 

Inclination 

P.A. 



Value 
Sc 

5.45' 

20 h 34 m 52 s .3 
60°09' 14.2" 
10 Mpc 
30° 
63° 



Reference 



de Vaucouleurs et al (1991) 



Tacconi & Young (1986) 
Corradi & Capaccioli (1991) 



Send offprint requests to: P.B.Alton, paul.alton@astro.cf.ac.uk 



Alton et al.: Submillimeter map of NGC 6946 



3 



In this paper, we address some of the issues outlined 
above by analysing data for NGC 6946. In most respects, 
this nearby (10 Mpc) galaxy resembles a typical Sc spi- 
ral and is oriented favourably for a study of face-on spiral 
structure (table Q). The centre of NGC 6946 evidently con- 
tains a large amount of recent star-formation (Tacconi & 
Young 1990; Engargiola 1991). However, there is no ev- 
idence from optical and near-infrared (NIR) imaging of 
tidal disturbance and indeed, the disk appears quite iso- 
lated on the sky (except for the low surface brightness 
spiral Cep 1, at 0.4 Mpc in projected distance, Burton 
et al. 1999). NGC 6946 is peculiar for its richness in CO 
emission, the derived molecular-to-atomic gas mass being 
quite high (Casoli et al, 1990). The density waves are not 
very regular nor contrasted, and in particular the total 
radio continuum emission is only slightly enhanced in spi- 
ral arms, indicating low compression (Boulanger & Vialle- 
fond, 1992). The contrast is even lower, when short spac- 
ings are added to interferometric maps (Frick et al 2001). 
Also, polarized radio emission is concentrated in interarm 
regions (Beck & Hoernes 1996). 

The global FIR-to-blue luminosity ratio is close to unity 
which is characteristic of quiescent spirals rather than an 
object undergoing extensive starburst activity (Soifer et 
al 1987; Bianchi et al 2000a). Our submm data for NGC 
6946 embrace a re-analysis of 850/im data collected by 
Bianchi et al (2000b) as well as some additional mapping 
observations, at the same wavelength, accessed from the 
SCUBA archive. We begin by examining the relationship 
between submm emission and neutral gas (HI and H2 ; §||) . 
A Monte Carlo radiative transfer model is then applied to 
a B-K colour map of NGC 6946 with a view to comparing 
visual optical depth with dust thermal emission 

(§§• The 

latter process leads to an upper limit in the grain emissiv- 
ity at 850/im (§||). Visual optical depth is then compared 
to column density of neutral gas in order to discover how 
the gas-to-dust ratio changes with galactocentric radius 
and gas-phase metallicity (§|j). Finally, we examine the 
distributions of HI, H2, submm thermal emission and vi- 
sual optical depth with respect to the spiral density wave, 
providing estimates of opacity and gas column density for 
both arm and interarm regions in NGC 6946 (§0). 

2. Observational Details 

Our submm data for NGC 6946 stem from two sources: 
22 hours of exposure taken by ourselves during April 
and June 1998, and a further 5 hours of observations 
retrieved from the SCUBA archive. The former have 
been presented briefly in Bianchi et al (2000b). All the 
data were reduced in a similar fashion to that described 
in Bianchi et al but more attention was now given to 
the subtraction of the sky signal (which is strong with 
respect to the source signal) . The Submillimeter Common 
User Bolometer Array (SCUBA) is mounted at the 
Nasmyth focus of the James Clerk Maxwell Telescope 



(JCMT) and provides simultaneous imaging at 450 and 
850/im by virtue of a dichroic beamsplitter (Holland et al 
1999), with a field of view of 2.3' . The shortwave array 
(450/im) consists of 91 bolometers (HPBW=7.5" ) whilst 
the longwave array (850/im) is composed of 37 elements 
(HPBW=14.7" ). At the shorter wavelength, we detect 
little more than the NGC 6946 nucleus and, therefore, 
our results will concentrate chiefly on the 850/mr emission. 

In scan-map mode, the telescope scans the source at 
a rate of 24 arcsec per second and at an angle allowing 
the source to be fully sampled (Pierce-Price et al 2000). 
The scan length was 9' , and spacing ~ 1' . During the 
scan, the secondary chops at a frequency of 7.8 Hz 
within the observed field facilitating the removal of sky 
background emission. Each scan-map is a convolution of 
the source with the chop and the profile of the source is 
eventually restored by a fourier transform (FT) technique. 
Scan-maps presented here cover most of the NGC 6946 
disk (-D25 — 10' ) and an individual scan took typically 
6 minutes. In order to sample emission over various 
scale-sizes, we adopted chop throws of 20, 30 and 65" 
along both the RA and Dec directions. 

The dedicated SCUBA software package, SURF 
(Jenness 1997), was employed to clean, flatfield and 
calibrate the images as well as correct for atmospheric 
attenuation. Since the chop takes place within the source 
field, the software also uses the edges of each scan to 
remove a linear baseline corresponding to the sky. At 
850/im, the zenith optical depth was usually between 
0.1 and 0.3. Scan-maps of Venus, processed in the 
same manner as the galaxy data, provided a calibra- 
tion from Volts to Janskys for each night's observations. 
The uncertainty in this calibration is estimated to be 15%. 

The Emerson II technique enabled a final image to 
be produced from the chop-convolved data (Emerson 
1995; Jenness et al 1998; Holland et al 1999). Scan-maps 
associated with each chop configuration were combined 
into a single image. Each of these 6 maps was then divided 
by the FT of the chop (a sine wave) in order to obtain the 
FT of the source. The FTs of the source were combined 
before the inverse FT process was carried out to produce 
the final image. Unfortunately, frequencies close to zero, 
corresponding to slow changes in the sky signal, are poorly 
sampled by this technique. This inadequacy is manifested 
by large-scale undulations in the convolved images for 
each chop configuration (Pierce-Price et al 2000). The 
difficulty was alleviated to some extent by adopting 
an iterative approach in the data processing. The final 
image at 850/im was fed back into an earlier stage of the 
reduction in order to help separate sky fluctuations from 
object signal. In our case, this was repeated 4 or 5 times 
before no further change was noted in the final image. 
Even after this process, some fluctuations of scalesize 
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~ 5' were still apparent in the final image. Attempts 
to remove these fluctuations by fitting a polynomial 
surface or using an unsharp-mask technique were largely 
successful and gave very similar results. We estimate 
that the largest detectable angular scale is of order 5' . A 
diffuse emission with a flat radial distribution larger than 
5' , the radius of NGC 6946, would not be detected. In 
the end, an unsharp-mask technique was selected. Here, 
the original image was smoothed with a gaussian of size 
5.5' FWHM and this smoothed version subtracted from 
the original version in order to give the signal above 
the "local" zero point. The resultant image was adopted 
as our final image at 850/xm. The la random noise in 
this final image is estimated to be 4.2 mJy/16"beam. 
The additional uncertainty associated with the sky 
subtraction (in particular the frequencies close to zero) 
is believed to be between 1.0 and 1.5 mJy/16" beam. 
This systematic error was estimated by dividing the 
scan-maps up into 3 sets (each set taken on different 
dates) , processing each data set separately and comparing 
the 3 resultant final images photometrically Similarly, 
various trials in the sky-subtraction process (polyno- 
mial surface, unsharp-masking etc.) suggest uncertainties 
of order 1 mJy/16" beam in establishing the zero baseline. 

For comparison with the SCUBA data, a 12 CO(2-l) 
image of NGC 6946 was kindly made available to us by 
Sauty et al (1998). These data were taken with the IRAM 
30-m telescope and have a spatial resolution of 13" . A 
new CO(l-O) map made with the 30-m telescope with 22" 
resolution will be soon available (Walsh et al, in press) 
and is discussed in Frick et al (2001). It is however very 
similar to that from Sauty et al (1998). Nicten et al. 
(1999) detected CO(4-3), which indicates the presence of 
warm gas, as far as the outer spiral arm. Bianchi et al 
(2000b) have already noted a strong correlation between 
the 850/zm dust thermal emission and the 12 CO(2-l) emis- 
sion line in NGC 6946. However, the relationship is not 
completely linear and there is also evidence for a bifurca- 
tion in scatter-plots of two quantities (discussed below). 
A very extensive image of NGC 6946 (26' x 26'), taken 
in the 21cm atomic hydrogen line, has been made avail- 
able to us more recently by F.Boulangcr. These data have 
a spatial resolution of 25" and combine both single-dish 
and interferometer data (Boulanger & Viallefond 1992). 
Deep images in the B and K wavebands, taken with the 
KPNO 0.9-m and 2.3-m WIRO telescopes, were obtained 
in April 1996 and November 1995 respectively (Trewhella 
1998). The K-band image was obtained by means of the 
128 x 128 pixel Michigan Infrared Camera and covers an 
area of 7.4' x 7.4' . 
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Fig. 1. SCUBA 850/im greyscale image (top) superim- 
posed by 850^m contours. Isophotes start at 3a and are 
separated by 0.5 mag. The beam (18" FWHM) is shown 
with white contours at the bottom-right of the plotted 
region. In the bottom figure, the submm isophotes have 
been projected onto a B-band image. 



3. Relationship between submm, CO and HI 
emissions 



The submm emission is a function of dust column density 
and dust temperature, while the HI and CO emissions 
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are function of gas column density and gas excitation. All 
three can be considered as gas tracers, with assumptions 
on gas-to-dust ratios, dust temperature and gas tempera- 
ture and density. 

The 850/im emission detected from NGC 6946 consti- 
tutes almost entirely thermal emission from dust grains 
(at T ~ 15K) heated by stellar radiation (Israel et al 
1999). This assumption begins to break down for the 
central 60" of the disk where contamination by the 
12 CO (3-2) line (A = 870/im) to the SCUBA longwave 
filter might be as high as 33% (Mauersberger et al 1999; 
Dumke et al. 2001). In fig. we show the SCUBA 850/im 
image after smoothing to a spatial resolution of 18" 
(FWHM). The salient features of this image are similar 
to those described in detail by Bianchi et al (2000b). The 
nuclear region is a strong source of submm radiation, with 
a maximum surface brightness of 250±40 mJy/18"beam 
at 850/^m. There are also local 'hot-spots' situated at 
~ 2.5' away from the nucleus and position angles of 44°, 
99° and 203°. These coincide with prominent HII regions 
in B-band and Ha images of NGC 6946. The inner disk 
(2.5' x 2.5' ) is reasonably well detected in our SCUBA 
image. 

To compare, on equal terms, the distributions of 
neutral gas and dust thermal emission, we smoothed 
the 12 CO (2-1) and 850/^m images to the same spatial 
resolution as the 21cm HI map (25" ) using a circular 
gaussian filter. Photometric sampling was then carried 
out, on all 3 images, for those regions where significant 
850/im emission (> 2a) had been detected. This was 
conducted using a 25" diameter aperture for a total of 
46 regions. To convert the 12 CO(2-l) line emission to 
molecular hydrogen column density, N(H 2 ), we adopted 
0.4 as the ratio of CO(2-l) to CO(1-0) intensity (as 
observed by Casoli et al 1990 in the main disk). We then 
assumed X = 1.5 x 10 20 cm~ 2 K kms -1 for the conversion 
of CO(1-0) to N(H 2 ). We recognize sizeable uncertainties 
in our conversion parameters, in particular a factor 2 
in the adopted 'X' factor. More insidious perhaps is the 
expected radial change in 'X' which is probably as large as 
any uncertainty in the global value (Neininger et al 1996; 
Dumke et al 1997). Furthermore, the ratio of CO(2-l) 
to CO(1-0) intensity, within the 25" of NGC 6946, has 
been shown to be closer to unity rather than the value 
of 0.4 adopted here (Weliachew et al 1988; Casoli et al 
1990). Fortunately, many of the inference we make in this 
paper are based on the main disk and exclude the central 
1' of the galaxy. It is this region which is likely to show 
dramatic changes in the conversion factors. Finally, we 
emphasize that we have assumed a fairly low value for X, 
consistent with recent gamma ray studies and the latest 
review articles of the CO(1-0) to N(H 2 ) conversion factor 
(Maloney 1990). 



Fig. |2| demonstrates how the submm emission relates 
to the hydrogen column density derived from CO and 
21cm HI maps. The top graph indicates a strong positive 
trend between molecular gas and 850/im surface bright- 
ness, F$ 50 . Fig. H also suggests little correspondance be- 
tween atomic gas and the submm emission we detect. For 
the disk as a whole this is intuitively correct because the 
850/im flux density is clearly concentrated towards the 
centre whereas NGC 6946, like so many other spirals, is 
known to possess a central depression in HI column den- 
sity (Boulanger & Viallefond 1992). Including the atomic 
gas in the total hydrogen column density does little to 
strengthen the correlation between neutral gas and -Fssoj 
primarily due to the dominance of N(_ff 2 ) at higher submm 
surface brightnesses. A least squares fit between total col- 
umn of neutral hydrogen, N(HI+H 2 ), and F 850 yields the 
following: 

N[H ^t H2) = (°- 115 ± °- 0052 ) x F «50 + (1.3 ± 0.4) (1) 
where N(HI+H 2 ) and F S50 are expressed in atoms/cm 2 
and mJy/25"beam, respectively. The two fitted coeffi- 
cients for the correlation with N(H 2 ) alone are 0.117 ± 
0.0057 for the slope and —0.5 ±0.4 for the zero point. The 
addition of HI slightly improves the correlation, but the 
effect is not very large, given the low HI column densities, 
with respect to the derived H 2 ones. In order to derive a 
gas-to-dust ratio from this relationship we are first com- 
pelled to relate i^so to the optical depth of dust in NGC 
6946. This step is carried out in the next section. 

4. Relationship between submm and reddening 

In this section, we relate thermal emission from dust at 
850/im to reddening evident in optical images of NGC 
6946. The deep B- and K-band images, alluded to in 
, were aligned and convolved to a common spatial 
resolution of 4.2" and the two images divided to produce 
a B-K colour image in magnitudes. In fig. we show this 
reddening image with 850//m contours superimposed. 
NGC 6946 is typical of a 'grand design' spiral galaxy with 
rather prominent spiral dust lanes (Trewhella 1998). We 
note the presence of extensive reddening ~ 2' north-west 
of the nucleus, while our 21cm HI and CO images give 
the impression of little atomic and molecular gas at this 
location. The filamentary structure in this part of the 
B-K map is more suggestive of reddening by dust rather 
than an absence of blue stars. 

Fig. H reveals that the alignment between reddening 
lanes and dust thermal emission is surprisingly poor. We 
have already identified some of the submm emission as 
hot-spots situated close to optically-bright HII regions. 
Thus some of the brightest regions at 850/im may well cor- 
respond to warm grains heated by young stars (T=35K; 
Bianchi et al 2000b). In contrast, optical extinction is 
not sensitive, in any direct sense, to grain temperature 
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although dust clouds situated in front of bright, blue star- 
forming regions will be detected quite readily in B-K 
maps. The situation contrasts strongly with the conclu- 
sions drawn for the nearby spiral NGC 891 (Alton et al 
2000). Here, radial profiles along the major axis of this 
galaxy revealed a fairly good correspondance between op- 
tical reddening and 850/im emission. Furthermore, there 
were strong indications from the thermal spectrum of 
NGC 891 that the submm is dominated by cold grains 
(15-20K) rather than dust surrounding HII regions (Alton 
et al 1998b). One possible explanation is that diffuse emis- 
sion is dominating in NGC 6946, and remains undetected 
because of the face-on orientation (while it is detected in 
the edge-on NGC 891). Radial profiles of NGC 6946, to 
be discussed in §[?], suggest a stronger link between visual 
optical depth and 850/im surface brightness than might 
be surmised from fig. |[ 

In addition to accentuating the absorption lanes, 
our B-K colour map highlights the presence of strong- 
reddening towards the centre of NGC 6946. This is 
accompanied by an increase in 850/im surface brightness; 
both properties, taken together, suggest an increase in 
dust concentration towards the galaxy nucleus. For the 
analysis that follows, we assume that the reddening in 
NGC 6946 is solely attributable to extinction by dust 
(and not, for example, caused by a colour change in the 
underlying stellar population). We adopt the observed 
B-K=3.1 colour at a radius of 4.7' (R ~ 0.9 x R 25 ) as the 
intrinsic, unreddened colour. This radius of 4.7' is chosen 
because extinction is small there, and it is the largest 
radius where colours can be estimated properly. We then 
assume that any increase in the colour index with radius 
(A(B — K)) arises from increasing optical depth. The 
results in radial variation (see fig. 8) do not a posteriori 
show the need for a more sophisticated approach. A 
Monte Carlo radiative transfer model is then employed 
to infer the visual optical depth, ry, from A(B — K). 

The assumption that we are measuring the true, 
underlying colour at a radius of R = 0.9 x i?25 is fairly 
robust - the optical depth is generally believed to be quite 
low near the R25 (ry ~ 0.2; Alton et al 2001). However, 
our assumption that an increase in A(B — K) is purely 
attributable to dust is probably incorrect, although this 
viewpoint has been proposed by many previous authors 
(Peletier et al 1995; Beckman et al 1996). Peletier et al 
(1995) find a statistical increase in the ratio of B-band 
radial scale-length to K-band radial scale- length for 
spiral disks more inclined to the line-of-sight. This effect 
is most readily explained in terms of dust reddening 
because the ratio would not change if the scale-lengths 
were determined by the stellar population gradient. In 
contrast, for a sample of 86 face-on spirals, de Jong (1996) 
concluded that the change in scale-length across several 
optical and NIR passbands could not be understood 
without gradients in star-formation and metallicity being 



considered as a major factor. The assumption in this 
paper that increased B-K colour originates purely from 
increased optical depth will lead us to overestimate 
the opacity of the NGC 6946 disk. Furthermore, since 
the stellar population along the spiral arms is usually 
considered younger, and therefore bluer, than adjacent 
intcrarm regions, we will tend to both overestimate 
the amount of interarm dust and, therefore, overlook 
the propensity for dust to clump along the spiral arms 
(Knapen & Beckman 1996). 

Fig. U provides a breakdown of the reddening 
measured in fig. [3] on a pixel- by-pixel basis. The 
azimuthally-averaged radial profile in B-K colour has 
been superimposed. We interpret A(B — K) by means of 
the Monte Carlo radiative transfer (RT) model developed 
by Bianchi et al (1996). This simulation assumes smooth 
distributions of both stars and dust which decline expo- 
nentially with both increasing galactocentric radius and 
height above the midplane. The light paths of optical 
and NIR photons are followed, through a succession of 
scattering interactions, until absorption takes place or the 
photon in question exits from the galaxy. By comparing 
the number of photons that escape from each part of 
the object with the quantity emitted by stars in the 
underlying disk, it is possible to relate observed colour to 
visual optical depth. 

For NGC 6946, we adopt a 'standard geometry' with 
dust grains extending to half the scale-height of the 
stellar disk and the dust radial scale- length fixed at 1.5 
times the stellar scale-length. Acceptable variations in the 
scale-height and scale-lengths of dust with respect to the 
stars (Trewhella 1997; Bianchi 1999) lead to uncertainties 
of ~ 33% in the implied visual optical depth over most 
of the disk (~ 50% for the central arcminute). An i? 1 / 4 
bulge was added to the stellar light so that, in total, bulge 
stars contributed ^ of the total optical+NIR luminosity 
of the galaxy. Finally, it should be noted that no account 
is taken of dust clumping in the model which normally 
would lead us to underestimate the true optical depth by 
50-100% (Witt & Gordon 2000; Bianchi et al 2000c). 

Rather than attempt a global fit between the observed 
and simulated B-K colours, each pixel in fig. ^ was treated 
independently. A series of simulations were run, with the 
face-on, visual optical depth through the centre of the 
galaxy model varying between 2 and 10. Overall, the 
technique worked well; a notable exception being within 
45" of the nucleus where the observed A(B — K) was too 
high to give a reliable indication of the underlying optical 
depth (B-K effectively saturates as ry approaches ~ 10). 

Using the method outlined above a map of visual 
optical depth was created to the same spatial resolution 
as our K-band image (4.2" FWHM). To compare with 
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the SCUBA image (next section), the Ty map was then 
smoothed to 18" , yielding a central, face-on opacity of 
t v — 7.6 and an exponential scale-length of 98" , close to 
that inferred from ISO 200/im images. We conclude that 
the visual optical depths derived from our RT model are 
reasonable. 

Fig. U shows how the column density of neutral hydro- 
gen in NGC 6946, N(HI + H 2 ), relates to both reddening, 
A(B — K), and the optical depth derived from RT mod- 
elling, Ty. If we exclude central regions saturated in B-K 
(and corresponding to Ty > 6), a strong linear relationship 
emerges as follows (correlation coefficient 0.91): 
N(HI + H 2 " 



10 21 



(1.0±0.1)tv + (2.1±0.3) 



(2) 



This is half the corresponding value measured within the 
solar neighborhood (N(HI+H 2 ) = 2.1 x 10 21 t v ; Bohlin et 
al 1978), and could be due to the difference in metallicity. 
To obtain a total dust mass, we now convert Ty into dust 
column density and integrate over the whole disk. To do 
this, we employ the relation given by Xilouris et al (1997) 
corresponding to 'classical' Milky Way grains of size 0.1/tm 
and material density 3 g cm -3 (Hildebrand 1983). Thus: 

M D = J 0.27t v (R)2ttR dR (3) 

where the integral is over radius R, in metres, along the 
disk. Here, t v (R) is the face-on optical depth at R and 
Md is the total dust mass in g. Integrating eq. [| out to 
a radius R ~ 4' , we derive M D = 8.6 x 1O 7 M for NGC 
6946. Integrating N(HI + H 2 ) over the same region leads 
us to a global gas-to-dust ratio of 90. This value is con- 
sistent with the somewhat shallower gradient derived in 
eq. ^| compared with the solar neighborhood. Given our 
fundamental assumption in the derivation of Ty (that B- 
K increases due to extinction only), our estimate of the 
gas-to-dust ratio must be considered a lower limit. 

Let us note that this value of 90 for the gas-to-dust 
ratio is about a factor 2 lower than the average of 200 
found towards 7 spiral galaxies (including NGC 6946) by 
Alton et al (1998a). This is due essentially to the lower 
X=H2/CO conversion ratio adopted here (X=1.5 instead 



of 2.8 10 2U cm~ 2 (K km s" 



) . Also the 850/xm SCUBA 



emission could trace even colder dust than the ISO 200/im 
emission used by Alton et al (1998a). For comparison, the 
gas-to-dust ratio for the solar neighborhood is 150, and 
100-300 for NGC 891, according to the adopted conversion 
ratio or dust emissivity (Guelin et al. 1993, Alton et al 
2000). 

5. Inferred submm emissivity 

The submm emission is not only a function of the dust 
column density and the dust temperature, but also of the 
nature of grains, which is reflected in their emissivity Q 
at various wavelengths. 

A comparison between the SCUBA image and the op- 
tical depth just derived from RT modelling should allow 



us to constrain the 850/im emissivity of dust grains in 
NGC 6946, Q(850/ito). The mass absorption coefficient 
(k), which is often used, is related to the grain emissiv- 
ity as follows: Ksso^m = 3 -^l|£^!Iii j where a and p are 
the radius and material density of the grains respectively 
(Hildebrand 1983). It can be shown (e.g. Alton et al 2000) 
that optical depth and submm brightness are related as 
follows: 
_ Q(V) 1.67 x 10~ 18 



cos{i) F 8 5o 



(4) 



v Q(850/im) 5(850/im, T) 
where t v is the face-on optical depth, Q(V) is the 
emissivity at A = 0.55/im, i is the inclination of the disk 
(i — 90° for edge-on), Fg5o denotes the surface brightness 
of the 850/im map, in Jy/18"beam, and B(850fim,T) 
refers to the blackbody intensity at 850/mi for grains 
of temperature T. We can see that an estimate of -fr— 

^850 



would serve as a handle on 



Q(V) 

Q(850m«j) ' 



With current observational limits there exists almost 
an order of magnitude uncertainty in Q(850/im) (and 
hence Q^^m) ' smce Q(Y) i s relatively well constrained 
to ~ 1.5; Whittet 1992). Indeed, emissivity values 
for wavelengths close to 1 mm are often extrapolated 
from observations of Galactic reflection nebulae near 
100-200/im (see Alton et al 2000 for a full discussion). 
At the same time we cannot overstate the importance 
of determining the submm/mm emissivity because, 
once Q{\) is known, submm/mm observations can then 
be furnished with visual optical depths in a relatively 
straightforward manner (Hildebrand 1983; Hughes et al 
1993). 

Radial profiles were fitted to both SCUBA and visual 
optical depth maps (both after smoothing to a spatial res- 
olution of 18" ). Values pertaining to the central arcminute 
of the object were discarded due to suspected saturation 
in the B-K colour. Fig. ^ shows that optical depth Ty 
compares rather well to submm surface brightness 2*850, 
when azimuthally averaged (while it was not true locally, 
cf fig. ||). A least-squares best fit between optical depth 
and submm brightness yielded the following: 

-p^ = 573 ±108 (5) 

Given the already mentioned uncertainties on Ty and 
2^850) this value of the ratio is an upper limit. If we are to 



derive 



Q(V) 



Q(850fj.m) fr 0111 ec L s - H an( i H we must have some no- 
tion of the grain temperature T. Unfortunately, no signal 
was recorded from the main disk of NGC 6946 during our 
450/im SCUBA observations (which would have allowed 
us to plot the thermal spectrum). ISO 200/im images of 
NGC 6946, in conjunction with IRAS lOO^rn HiRes data, 
suggest dust temperatures close to 20K in the main disk. 
RT modelling by Bianchi et al (2000a) points to T falling 
from 23K at the centre of NGC 6946 to 17K at the R 25 
radius (similar to values proposed for the diffuse Galactic 
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ISM; Mathis et al 1983; Draine & Lee 1984; Reach et 
al 1995; Masi et al 1995). Observations at both 450 and 
850/xm for the nearby spirals NGC 891 and NGC 7331 
suggest that 90% of galactic dust maintains a tempera- 
ture between 15 and 20K (Alton et al 1998b; Alton et 
al 2001; see also Guelin et al 1993). To translate ^ 

we adopt T = 18K. We note that a 50% 



into 



Q(V) 



Q(850jLim) 



Q(V) 



uncertainty in T will produce a 50% variation in Q(g 50A1Tra ) 
since i^so oc T in the Rayleigh- Jeans tail of the spec- 
trum. The substitution of temperature into eq. |] yields 
Q (850 (iL) — 1.4 x 10 5 and a mass absorption coefficient 
of K(850/im) ~ 0.26 cm 2 g~ 1 is we assume Q(V) = 1.5 
and grains of radius 0.1/im and material density 3 g cm~ 3 . 

The implied 850/im emissivity is a factor 3 lower than 
that measured by COBE for diffuse Milky Way dust (sug- 
gesting we detect too little 850/mi radiation compared 
with the optical depth inferred from RT modelling) . How- 
ever, we recall that our ry values are only upper limits as 
we assumed from the outset that B-K only changes due to 
extinction. The aforementioned 850/im emissivity should 
therefore be considered a lower limit. Also, there might 
be an underestimation of the emissivity, due to the pos- 
sible existence of spatially flat diffuse emission and the 
uncertainty in the background subtraction. In NGC 891, 
the submm emissivity was estimated 9 x higher than that 
for NGC 6946 (Alton et al 2000), which might be due to 
more thorough detection of the diffuse emission, because 
of edge-on orientation. 



6. Gas-to-dust Ratio 

Fig. U reveals how the derived gas content varies with re- 
spect to dust, with distance from the nucleus. 

In one case, the optical depth derived from RT 
modelling has been used to infer the dust mass per 
unit area (via eq. ||). The gas surface density includes 
both atomic and molecular gas (the latter based on 
X = 1.5 x 10 20 cuT 2 K kms" 1 ). The radial gas-to-dust 
ratio has also been calculated using the submm profile 
discussed in §|| and converting Fg^Q to ry with eq. pi 
The normalisation of -Fsso to Ty explains why the 
mean gas-to-dust ratio indicated by both gas-to-dust 
profiles in fig. ^ is more or less the same. Again, it must 
be emphasized here that both curves are only lower limits. 

We should remind here that the molecular gas is domi- 
nating in NGC 6946, and its amount is estimated from the 
CO emission, which is itself dependent on metallicity. In 
NGC 891 for example, both CO emission and submillime- 
ter dust emission have the same dependence with radius, 
suggesting that the CO emission falls also proportionally 
to the metallicity (Guelin et al. 1993). 

A dashed line in fig. || represents the gas-phase metal 
abundance as measured from emission-line spectra (Vila- 



Costas & Edmunds 1992). The vertical axes in the plot 
have been aligned so that solar metallicity corresponds to 
a gas-to-dust ratio of 150. The proximity of our gas-to- 
dust estimates to the dashed line is probably fortuitous 
given the factor 3-4 uncertainty in deriving the former (X 
factor, CO(2-1)/CO(1-0) etc.). 

Nearly no radial variation of the gas-to-dust ratio is 
observed. This means that the CO tracer is also dependent 
on metallicity, in a similar way than the dust tracer. 

7. Dust and gas within arm and interarm regions. 

The arm-interarm contrast in dust and gas is difficult 
to determine because of limited spatial resolution - at 
60/im, an IRAS beam of ~ 1' only begins to resolve spiral 
disks that are closer than ~ 5 Mpc. SCUBA, with a beam 
of 16" at 850/im, can probe spiral arms separated by 
1' (such as those in NGC 6946). However, the problem 
becomes one of sensitivity, with face-on disks possessing 
surface brightnesses of only a few mJy. 

Although our 21cm HI observations are somewhat too 
coarse for the task (25" FWHM), we have attempted 
a decomposition of the gas and dust associated with 
the spiral structure in NGC 6946, concentrating on the 
NE where the 850/mi emission is brightest. A series of 
8 cross-sections across the arm-interarm interface were 
produced, sampling the opacity as well as the submm, 
CO and 21cm HI emission. The transects were then 
stacked, assuming an average arm separation of 60" (A) 
and stretching or compressing to take into account the 
variation of A with galactocentric radius. 

Fig. H depicts the mean profiles generated in this way. 
The vertical dashed lines define arm and interarm regions 
on the basis of the K-band profile (believed to trace the 
underlying stellar mass density). All the tracers of dust 
and gas plotted in fig. ^ show some modulation with 
the spiral density wave except perhaps the optical depth 
derived from RT modelling. The lack of variation here 
probably stems from the assumption that an identical 
intrinsic B-K colour exists for the underlying stellar 
population in both arm and interarm regions. In reality, 
the arms will tend to be inherently bluer on average, due 
to the presence of young stars, leading us underestimate 
the arm opacity with respect to the interarm value. 
The two forms of neutral hydrogen seem complementary 
in the sense that a reduced presence of HI along one 
arm is balanced by an enhanced level of H2 and vice versa. 

Table || contains estimates of the optical depth and gas 
column density based on the profiles presented in fig. ^[ 
The visual opacity is calculated variously from: radiative 
transfer (row 3); 850/im surface brightness using eq. |^ (row 
13); and the gas column density N(HI + H2) (rowl2), 
the last of these assuming tb — Ty — R ^tt^^t-^H — 7 
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Table 2. Properties of the arm and interarm regions in 
NGC 6946. Measurements are shown for the following: 
K-band light (row 1), 850/xm surface brightness (row 2), 
visual optical depth (row 3), 12 CO(2 — 1) intensity (row 
4), H2 column density (row 5), intensity in the 21cm HI 
line (row 6), HI column density (row 7), total hydrogen 
column density (row 8), ratio of atomic to molecular col- 
umn density (row 9), visual optical depth inferred from 
H2 (row 10), visual optical depth inferred from HI (row 
11), visual optical depth inferred from H2+HI (row 12) 
and visual optical depth inferred from the 850/im sur- 
face brightness (row 13). The H2 column density assumes 
X = 1.5 x 10 20 cm~ 2 K kms -1 and 0.4 for the ratio of 
CO(2-l) to CO(l-O) intensity. For inferences of visual op- 
tical depth see text. Due to the cancellation of system- 
atic errors, uncertainties in the interarm-to-arm ratios are 
sometimes proportionally smaller than those pertaining to 
the arm or interarm region alone. Note that none of the 
quantities appearing in the table have been corrected to 
face-on values. 



tively, for the optical depth of the interarm and arm re- 
gion at a radius of 2.3' (table ||). In addition, we assume 
a radial scale-length of 98" for both arm and interstitial 
environments (§[!]). We suppose that interarm regions oc- 
cupy only 50% of the disk at R = 0.5 x R25, increasing 
linearly with R from 0% at the nucleus to 100% at the 
i?25 radius. Integrating eq. || suggests that at least 15% of 
galactic dust is located between the spiral arms. A similar 
calculation for HI+H2 (row 8, table |J) indicates that at 
least 68% of the neutral gas in NGC 6946 can be consid- 
ered as lying outside the spiral arms. 

8. Summary 

We have analysed 850/im SCUBA images of the nearby, 
(almost) face-on spiral galaxy NGC 6946. The submillime- 
ter brightness level is low over most of the disk (often a 
few mjy at 850/im) and, therefore, care has to be exer- 
cised in the subtraction of sky emission. Our main results 
can be summarised as follows: 



Quantity Units 



K-band 

-£850 
TV 

Ico 
N(H 2 ) 

1 21cm 

N(HI) 
7V(i/ 2 ±HI) 

N(H 2 ) 
N(HI) 

tv{H 2 ) 
TV {HI) 
7V(#2+HI) 

iv(Fgso) 



counts 
mJy/18" 

K km/s 
10 21 cm- 2 
K km/s 
10 21 cm- 2 
10 21 cm- 2 



Arm 

average 

13.4 ± 1.2 
6.4 ±1.5 
< 2.05 
5.32 ± 0.4 
2.00 ±0.1 
941 ± 58 
1.79 ±0.1 
3.8 ±0.2 
1.12 ±0.1 
0.96 ± .08 
0.86 ± .04 
1.8 ±0.1 
3.7 ±0.9 



Interarm 
average 

11.4 ± 1.2 
2.7 ± 1.5 
< 2.04 
4.10 ±0.4 
1.54 ±0.1 
787 ± 58 
1.50 ±0.1 
3.0 ±0.2 
1.03 ±0.1 
0.74 ± .08 
0.72 ± .04 
1.5 ±0.1 
1.5 ±0.9 



Interarm 
/ Arm 

0.86 ± .02 
0.42 ±0.3 
< 1.0 
0.77 ± .07 
as above 
0.84 ± .05 
as above 
0.78 ± .07 



1. 



as observed in the solar neighborhood (Bohlin et al 
1978). Our estimate of ry based on RT modelling is 
expected to be an overestimate in the interarm region. 
Surprisingly, however, the gas column in such interstitial 
regions suggests equally high values of Ty. The submm 
observations favour quite a high interarm opacity as well 
(1.5 ± 0.9) but the combination of low signal-to-noise 
and difficulties in relating grain thermal emission to dust 
column density make the uncertainties here considerable. 
Attempts to determine interarm opacity by analysing 
the colour, surface brightness or number of background 
galaxies visible through a foreground spiral disk have 
generally led to values of less than unity for the face-on 
optical depth (White & Keel 1992; White et al 1996; 
Berlind et al 1997; Gonzalez et al 1998). 

Although the definition of arms and interarms regions 
is subjective, we have attempted to compute the gas and 
dust mass in both regions. We adopt 1.5 and 3.7, respec- 



Comparing the 850/im flux with 12 CO(2-l) and 21cm 
HI line intensity indicates a good correlation between 
dust thermal emission and CO emission. No correlation 
was detected between submm emission and atomic gas 
emission. 

We place an upper limit on the visual optical depth 
in NGC 6946 (ry) by means of a B-K colour map of 
the disk. We assume that the B-K colour at a radius 
0.9i?25 corresponds to the underlying intrinsic colour 
of the stellar disk, in both arm and interarm regions, 
and that the observed increase in B-K with diminishing 
radius is attributable solely to extinction by dust. With 
these conditions we carry out a Monte Carlo radiative 
transfer simulation of the disk, inferring visual optical 
depth from B-K colour. 

The resultant map in Ty relates fairly well to the distri- 
bution of neutral gas (HI±H2). There is no significant 
radial variation of the derived gas-to-dust ratio, which 
is easy to understand, given that the estimated gas 
content is dominated by the molecular gas, therefore 
is estimated through the CO emission tracer. Since the 
CO emission is also dependent on metallicity, it is not 
possible to derive the true gas-to-dust ratio. This will 
await a large-scale mapping of other tracers, such as 
the pure rotational lines of molecular hydrogen (e.g. 
Valentijn & van der Werf, 1999). 

Radial profiles of our optical depth map and SCUBA 
image indicate an association between Ty and (dif- 
fuse) submm emission, F$5q. A linear regression yields 
-pt— = 573 ± 110 and, by corollary, a 850/im emissivity 
which is no more than 3 times lower than the COBE 
value for Galactic grains (Boulanger et al 1996). The 
emissivity is also 9 times lower than in the edge-on 
galaxy NGC 891. This is however a lower limit, given 
the uncertainties on the visual opacity, and the possi- 
bility of underestimating the large-scale submm emis- 
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sion. The latter could be estimated in the future by 
more sensitive and more extended submm arrays. 
5. A decomposition of spiral structure in NGC 6946 sug- 
gests an interarm optical depth of between 1 and 2 at a 
mean radius of 0.4i?25- These surprisingly high values 
represent 40-80% of the visual opacity along the arms. 
Furthermore, by a similar analysis, about two-thirds 
of the neutral gas could be considered as lying outside 
the spiral arms. 

Acknowledgements. It is a pleasure to thank Gerald Moriarty- 
Schieven and Tim Jenness for their support during observa- 
tions and data reduction. We are grateful to Steve Eales and 
Loretta Dunne for stimulating discussions on grain opacity and 
related matters, and to the referee for his careful reading and 
positive criticism. 

References 

Agladze, N., Sievers, A., Jones, S., Burlitch, J., Beckwith, S., 

1994, Nature, 372, 243 
Alton, P.B., Trewhella, M., Davies, J., 1998a, A&A, 335, 807 
Alton, P., Bianchi, S., Rand, R., Xilouris, E., Davies, J., 

Trewhella, M., 1998b, ApJL, 507, L125 
Alton, P., Xilouris, E., Bianchi, S., Davies, J., Kylafis, N., 2000, 

A&A, 356, 795 
Alton, P.B., Lequeux, J., Bianchi, S., 2001, A&A, 366, 451 
Beck R., Hoernes P.: 1996, Nature 379, 47 
Beckman, J., Peletier, R., Knapen, J., Corradi, R., Gentet, L, 

1996, ApJ, 467, 175 
Berlind, A., Quillen, A., Pogge, R., Sellgren, K., 1997, AJ, 114, 

107 

Bianchi, S., Ferrara, A., Giovanardi, C, 1996, ApJ, 465, 127 
Bianchi, S., 1999, PhD Thesis, University of Cardiff 
Bianchi, S., Davies, J., Alton, P., 1999, A&A, 344, LI 
Bianchi, S., Davies, J., Alton, P., 2000a, A&A, 359, 65 
Bianchi, S., Davies, J., Alton, P., Gerin, M., Casoli, F., 2000b, 
A&A, 353, L13 

Bianchi, S., Ferrara, A., Davies, J., Alton, P., 2000c, MNRAS, 

311, 601 

Block, D., Witt, A., Grosbol, P., Stockton, A., Moneti, A., 

1994, A&A, 288, 383 
Bohlin, R., Savage, B., Drake, J., 1978, ApJ, 224, 132 
Boulanger, F., Viallefond, F., 1992, A&A, 266, 37 
Boulanger, F., Abcrgcl, A., Bernard, J. P., et al, 1996, A&A, 

312, 256 

Burton, W. B., Braun, R., Walterbos, R. A. M., Hoopes, C. 

G.: 1999, AJ 117, 194 
Casey,S., 1991, ApJ, 371, 183 

Casoli, F., Clausset, F., Combes, F., Viallefond, F., Boulanger, 

F., 1990, A&A, 233, 357 
Corradi, R., Beckman, J., Simonneau, E., 1996, MNRAS, 282, 

1005 

de Jong, R., 1996, A&A, 313, 377 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H., et al, 1991, 
Third Reference Catalogue of Bright Galaxies (New York, 
Springer) 

Draine, B., Lee, H., 1984, ApJ, 285, 89 

Dunne, L., Eales, S., Edmunds, M., et al, 2000, MNRAS, 315, 
115 

Dumke, M., Braine, J., Krause, M., et al, 1997, A&A, 325, 124 



Dumke, M., Nieten, Ch., Thuma, G., Wielebinski, R., Walsh, 
W.: 2001, A&A, 373, 853 

Edmunds, M., Eales, S., 1998, MNRAS, 299, L29 

Emerson, D., 1995, in Astronomical Society of Pacific Confer- 
ence Series, Vol. 75 

Engargiola, G., 1991, ApJS, 76, 875 

Fall, M., Pei, Y., 1993, ApJ, 402, 479 

Frick P., Beck R., Berkhuijsen E.M., Patrickeyev I.: 2001, MN- 
RAS 327, 1145 
Genzel, R., Cesarsky, C.J.: 2000, ARAA 38, 761 
Gonzalez, R, Allen, R., Dirsch, B., et al, 1998, ApJ, 506, 152 
Guelin, M., Zylka, R., Mezger, P., et al, 1993, A&A., 279, L37 
Guelin, M., Zylka, R., Mezger, P., Haslam, C, Kreysa, E., 1995, 

A&A, 298, L29 
Heisler, J., Ostriker J., 1988, ApJ, 332, 543 
Hildebrand, 1983, QJRAS, 24, 267 

Holland, W., Robson, E., Gear, W. et al, 1999, MNRAS, 303, 
659 

Hughes, D., Robson, E., Dunlop, J., Gear, W., 1993, MNRAS, 
263, 607 

Israel, F., van der Werf, P., Tilanus, R., 1999, A&A, 344, L83 
Jenness, T., 1997, Starlink User Note 216.1 
Jenness, T., Lightfoot, J., Holland, W., 1998, Proc. SPIE, 3357, 
548 

Knapen, J., Beckman, J., 1996, MNRAS, 283, 251 

Maloney 1990, in Interstellar Medium of Galaxies ed. Thron- 

son, H., Shull, J. 
Masci, F., Webster, R., 1999, MNRAS, 305, 937 
Masi, S., Aquilini, E., Boscaleri, A., et al, 1995, ApJ, 452, 253 
Mathis, J., Mezger, P., Panagia, N., 1983, A&A, 128, 212 
Mauersberger R., Henkel, C, Walsh, W., Schulz, A., 1999, 

A&A, 341, 256 

Neininger, N., Guelin, M., Garcia-Burillo, S, et al, 1996, A&A, 
310, 725 

Nieten, Ch., Dumke, M., Beck, R., Wielebinski, R.: 1999, A&A 
347, L5 

Peletier, R., Valentijn, E., Moorwood, A., et al, 1995, A&A, 
300, LI 

Pierce-Price, D., Richer, J., Greaves, J., et al, 2000, ApJ, 545, 
L121 

Rand, R., Lord, S., Higdon, J., 1999, ApJ, 513, 720 

Reach, W., Dwek, E., Fixsen, D., et al, 1995, ApJ, 451, 188 

Rengarajan, T., 1984, A&A, 140, 213 

Sauty, S., Gerin, M, Casoli, F., 1998, A&A, 339, 19 

Seigar, M., James, P., 1998, MNRAS, 299, 685 

Siebenmorgen, R., Kruegel, E., Chini, R., 1999, A&A, 351, 495 

Soifer, B., Sanders, D., Madore, B., et al, 1987, ApJ, 320, 238 

Sopka, R., Hildebrand, R., Jaffe, D., et al, 1985, ApJ, 294, 242 

Tacconi, L., Young, J., 1986, ApJ, 308, 600 

Tacconi, L., Young, J., 1990, ApJ, 352, 595 

Trewhella, M., 1997, PhD Thesis, University of Cardiff 

Trewhella, M., 1998, MNRAS, 297, 807 

Valentijn E.A., van der Werf P.: 1999, ApJ 522, L29 

Vila-Costas, M., Edmunds, M., 1992, MNRAS, 259, 121 

Weliachew, L, Casoli, F., Combes, F., 1988, A&A, 199, 29 

White, R., Keel, W., 1992, Nature, 359, 129 

White, R., Keel, W., Conselice, C, 1996, in New Extragalac- 

tic Perspectives in the New South Africa, edt. Block, D., 

Greenberg, M. (Kluwer) 
Whittet, D., 1992, Dust in the Galactic Environment (IOP 

Publishing) 



Alton et al.: Submillimeter map of NGC 6946 



11 



Williams, D., Taylor, S., 1996, QJRAS, 37, 565 
Witt, A, Gordon, K., 2000, ApJ, 528, 799 
Xilouris, E., Kylafis, N., Papamastorakis, J., Paleologou, E., 
Haerendel, G., 1997, A&A, 325, 135 




100 200 300 

F 85o (mJy/25"beam) 




100 200 300 

F 85o (mJy/25"beam) 




100 200 300 

F 850 (mJy/25"beam) 



Fig. 2. Relationship between neutral gas and 850^<m emis- 
sion in NGC 6946. The submm surface brightness is de- 
noted by F 850 whilst N(HI) and N(H 2 ) refer to the col- 
umn density of atomic and molecular hydrogen respec- 
tively. In all cases, an error-bar at the top-left of the plot- 
ted region indicates the random error. The inset in the top 
graph magnifies data points occurring at lower 850/im sur- 
face brightness. 
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Fig. 3. B-K colour image of NGC 6946 displayed in 
greyscale (reddened regions appear as white). This is quite 
similar to the extinction map of Trewhella (1998) combin- 
ing optical and far-infrared. The 850/im contours appear- 
ing in fig. |l| have been superimposed. 
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Fig. 4. A breakdown of the B-K colour in fig. [| Dots 
represent the difference in colour of each pixel, A(B — K) 
with respect to the colour at the disk edge (R = 4.7'). 
The baseline at R — 4.7' has been established by means 
of an azimuthally- averaged profile (solid circles and solid 
line). 




0.5 1 1.5 
A(B-K) (mag) 



Fig. 5. Neutral hydrogen column density (H2+HI) plot- 
ted against both B-K colour and visual optical depth, ry . 
The latter has been inferred from a Monte Carlo radia- 
tive transfer simulation. A least-squares linear fit between 
gas column density and ry is represented by the solid 
line (excluding ry > 7). The corresponding relationship 
in the Milky Way is also linear but with a gradient ap- 
proximately twice that shown for NGC 6946. 
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F 85o (mJy/18"beam) 

Fig. 6. The relationship between visual optical depth (ry ) 
and 850 ^m surface brightness (F 850 ) for NGC 6946. Er- 
rorbars are only shown for i<850 as the values for ry are 
almost certainly upper limits. The solid line shows a least- 
squares linear fit to the data (see text). 
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Fig. 7. The emissivity of dust grains in the FIR/submm 
waveband. Measurements are plotted for Galactic re- 
flection nebulae (Hildebrand 1983; Casey 1991), deeply- 
embedded FIR sources (Rengarajan 1984), dust envelopes 
of evolved stars (Sopka et al 1985), the edge-on spiral NGC 
891 (Alton et al 2000) and NGC 6946 (upper limit, cor- 
responding to our lower limit on submillimeter emissivity 
given in this paper). Agladze et al (1994) refers to lab- 
oratory experiments conducted on silicate particulates. 
Draine & Lee (1984) corresponds to a model for diffuse 
Galactic dust partly based on empirical data. We also 
show values derived from COBE and IRAS data based 
on emission from high-lattitudc Galactic dust (Boulangcr 
et al 1996 and Bianchi et al 1999 respectively). The dot- 
ted and dashed lines denote, respectively, Q(X) oc A -1,5 
and Q(A) oc A -2 and are arbitrarily positioned to pass 
through the Draine & Lee 125^m point. The right hand 
axis shows the equivalent mass absorption coefficient (k) 
assuming a radius and material density of 0.1/im and 3 
gcm~ 3 , respectively, for the dust grains and Q(V)~1.5 
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Fig. 8. The radial behaviour in gas-to-dust ratio (left axis) 
with the dust mass derived from ; (i) RT modelling (solid 
line) and (ii) 850/im emission (dotted line). The dashed 
line shows the gas-phase metallicity in NGC 6946 and its 
associated error (right axis, plotted as ^§-)- The vertical 
axes have been aligned so that a gas-to-dust ratio of 150 
coincides with solar metallicity. A vertical error-bar shows 
the probable systematic error in the metallicity measure- 
ment (Vila-Costas & Edmunds 1992). 




Fig. 9. Profiles of gas and dust across the arm and in- 
terarm regions of NGC 6946. In all cases the horizontal 
axis denotes the distance from the interarm centroid (r) 
in terms of the mean arm-to-arm separation (A). Dashed 
vertical lines delineate both the arm regions (r = — 0.75A 
tor = -0.25A and r = +0.25A to r = +0.75A) and the 
interarm region (r = — 0.25A to r — +0.25A) as defined 
on the basis of the K-band image (top- left). Profiles are 
shown for 850/um surface brightness {F 850 ), visual optical 
depth (ry), 12 CO(2-l) intensity (Ico) an d intensity in the 
21cm HI line (/21cm)- Errorbars for these quantities are 
only shown when the uncertainties are much larger than 
the plotted markers. The spatial resolution of the data is 
represented by a horizontal error-bar at the bottom-left of 
each panel. 



